Diffusive and localization behavior of electromagnetic waves in a two-dimensional 

random medium 
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In this paper, we discuss the transport phenomena of electromagnetic waves in a two-dimensional 
random system which is composed of arrays of electrical dipoles, following the model presented 
earlier by Erdogan, et al. (J. Opt. Soc. Am. B 10, 391 (1993)). A set of self-consistent equations 
is presented, accounting for the multiple scattering in the system, and is then solved numerically. A 
strong localization regime is discovered in the frequency domain. The transport properties within, 
near the edge of and nearly outside the localization regime are investigated for different parameters 
such as filling factor and system size. The results show that within the localization regime, waves 
are trapped near the transmitting source. Meanwhile, the diffusive waves follow an intuitive but 
expected picture. That is, they increase with travelling path as more and more random scattering 
incurs, followed by a saturation, then start to decay exponentially when the travelling path is large 
enough, signifying the localization effect. For the cases that the frequencies are near the boundary 
of or outside the localization regime, the results of diffusive waves are compared with the diffusion 
approximation, showing less encouraging agreement as in other systems (Asatryan, et al., Phys. 
Rev. E 67, 036605 (2003).) 

PACS numbers: 42.25.Hz, 41.90.+e 
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I. INTRODUCTION 

Waves are a ubiquitous phenomenon that is relevant 
to our everyday's life. Our knowledge about the Na- 
ture is mainly obtained through either acoustic or elec- 
tromagnetic waves. When waves propagate through a 
random medium, some peculiar properties would emerge 
and some of them have stood as a long outstanding prob- 
lem for physicists. The most intriguing and still un- 
solved puzzle perhaps is the so called Anderson localiza- 
tion. The concept of Anderson localization was initially 
proposed nearly half a century ago for the possible phe- 
nomenon of disorder-induced metal-insulator transition 
in electronic systems pj. To make it simple, Anderson 
localization refers to situations that when released in a 
random medium which could be, for instance, a free space 
with random potentials, electrons may go nowhere but 
stay close to the initial place. The envelop of the elec- 
tronic wave function subsequently reveals as an exponen- 
tial decay along any direction from the emitting point ; 
the length measuring the e-fold decay is called the local- 
ization length. The mechanism behind this property has 
been attributed purely to sufficient multiple scattering of 
electrons by the random potentials, a feature of the wave 
nature of electrons. 

Since its inception, the localization concept has opened 
a wide door for scientists from various backgrounds, and 
stimulated a vast body of research. The fact that elec- 
tronic localization is due to the wave nature of electrons 
is particularly important, as the concept may also be 
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applied to classical wave systems by analogy ;3]. The 
concept of localization has far reached many fields such 
as seismology^, oceanologyjfjj, and random lasers |g, to 
name just a few. The great efforts have been summarized 
in a number of reviews (e. g. |1 EE IS II H 111 111 03 ) 

Over the past two decades, localization of electro- 
magnetic and acoustic waves has been and continued to 
be a particularly attracting proble m, leading to a great 
amount of publicatio ns ( e. 
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sical waves has been detailed in the textbook^]. In spite 
of the tremendous efforts, however, there are still some 
open questions for wave localization in two dimensional 
(2D) random systems. For example, one of the questions 
concern with the transport behavior of waves within and 
outside localization regimes. It has been discussed much 
in the literature that within the localization regime the 
wave intensity follows a diffusive behavior for transport 
path smaller than localization length, whereas it follows 
an exponential decay when the path is greater than the 
localization length |^,|29j. In addition, it has been con- 
jectured that all waves are bound to be localized in 2D 
for any given amount of disorders or randomness 30], a 
conclusion significantly influencing nearly all later inves- 
tigations. If this general hypothesis is true, one would 
expect that the wave intensity would have to follow the 
diffusion to localization transition as the transport dis- 
tance increases although such a transition has not been 
solidly confirmed yet. Recent studies, however, tended 
to s ugg est that this conjecture may be generally invalid 
[in],|32j- One one hand, in 31], it has been suggested 
that acoustic waves could not be localized in certain fre- 
quency ranges. One the other, the research reported in 
[32( shows different frequency regimes of electromagnetic 
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wave transport. In certain ranges, the wave transport be- 
havior can be entirely described by the usual Boltzmann 
diffusion approximate theory, while in other ranges the 
behavior does not reveal a diffusive feature. The au- 
thors called the latter phenomenon as anomalous diffu- 
sion, and inferred as the incipience of Anderson localiza- 
tion. These works seem to indicate that as the frequency 
varies, there will be a transition between localization and 
non-localization even in two dimensions. Whether these 
observations also hold in a general perspective remains 
to be known. 

The difficulties in the study of possible non-localization 
or localization of EM waves mainly lie in a number 
of problems. First, wave localization only appears for 
strongly scattering media, and such a medium is often 
hard to find. Second, localization effects are often entan- 
gled with other effects such as dissipation, wave deflec- 
tion, or boundary effects[l3<], making data interpretation 
often ambiguous. Third, many two dimensional systems 
are exactly solvable by numerical computation. But the 
numerical simulation is very time-consuming and is obvi- 
ously limited by computing facilities with regard to some 
unavoidable problems such as finite size. How to find 
a suitable model that could ease these concerns poses a 
challenge problem in its own right. 

In a recent communication, a simple but seemingly re- 
alistic model system has been proposed to study EM lo- 
calization in 2D random media[33]. This model origi- 
nated from the previous study of the radiative effects of 
the electric dipoles embedded in structured cavities |3J|. 
It was shown that EM localization is possible in such a 
disordered system. When localization occurs, a coherent 
behavior appears and is revealed as a unique property 
differentiating localization from either the residual ab- 
sorption or the attenuation effects. 

With the present paper, we wish to explore further the 
transport properties of the system outlined in [33[ ■ The 
advantage of this system not only relies on its simplicity, 
but on its less time consumption. It has been shown that 
there is a single strong localization region in the system. 
We will investigate the transport behavior of wave inten- 
sity within, near or outside the strong localization regime. 
We study how the transport behavior depends on scat- 
terer's filling parameter, frequency, and sample size and 
so on. When the localization effect is less obviously, we 
compare the numerically evaluated wave intensity with 
the result from the the diffusion theory. 

The paper is structure as follows. The description of 
the system and the relevant theoretical modelling are pre- 
sented in next section, followed by the numerical sim- 
ulation and detailed discussion of results. A summary 
concludes the paper in the last section. 



II. THE SYSTEM AND THEORETICAL 
FORMULATION 

Here we present the system and the theoretical formu- 
lation. Although these have already been presented be- 
fore [3{|, for the sake of the convenience on the reader's 
part and easy discussion on our part, we purposely repeat 
the essence here. 



A. The system 

Following Erdogan et al. |34j . we consider 2D dipoles 
as an ensemble of harmonically bound charge elements. 
In this way, each 2D dipole is regarded as a single dipole 
line, characterized by the charge and dipole moment per 
unit length. Assume that N parallel dipole lines, aligned 
along the z-axis, are embedded in a uniform dielectric 
medium and randomly located at = 1, 2, ... , N). The 
cross section of the dipoles are assumed to be in the x — y 
plane. The averaged distance between dipoles is d. A 
stimulating dipole line source is located at r Sl transmit- 
ting a continuous wave of angular frequency u. By the 
geometrical symmetry of the system, we only need to 
consider the z component of the electrical waves. 

B. The formulation 

Although much of the following materials can be re- 
ferred to in [33|. we repeat the important parts here for 
the sake of convenience and completeness. 

Upon stimulation, each dipole will radiate EM waves. 
The radiated waves will then repeatedly interact with the 
dipoles, forming a process of multiple scattering. The 
equation of motion for the z-th dipole is 

d 2 q 2 d 

fori = 1,2,..., AT. (1) 

where wo, t is the resonance (natural) frequency, pi, qi 
and rrii the dipole moment, charge and effective mass per 
unit length of the i-th dipole respectively. E z {fi) is the 
total electrical field acting on dipole pi, which includes 
the radiated field from other dipoles and also the directly 
field from the source. The factor &o,i denotes the damping 
due to energy loss and radiation, and can be determined 
by energy conservation. Without energy loss (to such as 
heat), &o,i can be determined from the balance between 
the radiative and vibrational energies, and is given as|34| 



with e being the constant permittivity and c the speed of 
light in the medium separately. 

Equation Q is virtually the same as Eq. (1) in [34| . 
The only difference is that in [34| , E z is the reflected field 
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at the dipole due to the presence of reflecting surround- 
ing structures, while in the present case the field is from 
the stimulating source and the radiation from all other 
dipoles. 

The transmitted electrical field from the continuous 
line source is determined by the Maxwell equations [3^ 



if 



c 2 dt 2 



G (r-r s ) = -4^ w 2 po7T<5 (2) (r-r s )e~ 



where u> is the radiation frequency, and po is the source 
strength and is set to be unit. The solution of Eq. J3J is 
clearly 



G (f- f* a ) = {pou 2 po)i'fH^'\k\f- r a \)e~ 



(4) 



with k = iv/c, and Hq being the zero-th order Hankel 
function of the first kind. 

Similarly, the radiated field from the i-th dipole is 
given by 



<:) 2 



c 2 dt 2 



Gi{r- n) = Mo ^«5(2)(f - f-). (5) 



The field arriving at the i-th dipole is composed of the 
direct field from the source and the radiation from all 
other dipoles, and thus is given as 



N 



E z (r t ) = G (n - r s ) + £ G^n-r,). (6) 
Substituting Eqs. (3J|, ©, and J§| into Eq. QJ, and 



writing p l = p^e' 



, we arrive at 

(-uj 2 + loI^ - iujb 0ii )pi 



rrii 



N 



G (r. l -r s ) + J2 



iH^\k\fi -fj\)pj 



•CO 



This set of linear equations can be solved numerically for 
Pi. After pi are obtained, the total field at any space 
point can be readily calculated from 



N 



E z (r) = G (r-r s ) + J2G J (r 



rj). 



(8) 



In the standard approach to wave localization, waves 
are said to be localized when the square modulus of the 
field |-E(r)| 2 , representing the wave energy, is spatially 
localized after the trivial cylindrically spreading effect is 
eliminated. Obviously, this is equivalent to say that the 
further away is the dipole from the source, the smaller its 
oscillation amplitude, expected to follow an exponentially 
decreasing pattern. 

An alternative two dimensional dipole model was de- 
vised previously . The authors derived a set of linear 
algebraic equations, which is similar in form to the above 
Eq. JfJJl. However, the fundamental difference between 



the two models is with respect to the relation between 
the incident and scattered waves. In , the interaction 
between dipoles and the external field is derived by the 
energy conservation, while in the present case the cou- 
pling is determined without ambiguity by the Newton's 
second law. The former leads to an undetermined phase 
factor. According to, e. g. Refs. 0,|3(|, the energy con- 
servation can only give the radiation factor in Eq. @. 

There are several ways to introduce randomness to 
Eq. J7J. For example, the disorder may be introduced by 
randomly varying such properties of individual dipoles as 
the charge, the mass or the two combined. This is the 
most common way that the disorder is introduced into 
the tight-binding model for electronic waves[37|. In the 
present study, the disorder is brought in by the random 
distribution of the dipoles. 

For simplicity yet without losing generality, assume 
that all the dipoles arc identical and they are randomly 
distributed within a square area. The source is located at 
the center (set to be the origin) of this area. For conve- 
nience, we make Eq. Q non-dimensional by scaling the 
frequency by the resonance frequency of a single dipole 
luq. This will lead to two independent non-dimensional 

parameters b — and b' Q = ^ (u^J ' B°th parame- 
ters may be adjusted in the experiment. For example, the 
factor bo can be modified by coating layered structures 
around the dipoles[23|. Then Eq. Q becomes simply 



(-/ 2 + 1 - ib' Q ) Pi = 



ibf 



N 



p H^\k\f i -f s \)+ J2 PjH^m-rjl) 

3 = ^,j¥=i 



(9) 



with / = -jj-. Eq. © is self-consistent and can be solved 
numerically for pi and then the total field is obtained 
through Eqs. ©, © and JSJ|. 



C. Expected transport behavior 

1. In general 

Following [38|, a general consideration of the spatial 
behavior of wave transport is possible. Consider a wave 
transmitted in a random medium. The transport equa- 
tion for the total energy intensity /, i. e. < \E\ 2 >, may 
be intuitively written as 



dl 

— = -al, 
ax 



(10) 



where a represents decay along the path traversed. After 
penetrating into the random medium, the wave will be 
scattered by random inhomogeneities. As a result, the 
wave coherence starts to decrease, yielding the way to 
incoherence. Extinction of the coherent intensity Ic, i- e. 
< E > | 2 , is described by 



dl c 
dx 



-lie, 



(11) 
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with the attenuation constant 7. Eqs. (|10|) and (|llf> lead 
to the exponential solutions 

I(x) = l(0)e' ax , and I c {x) = I{Q)e^ x . (12) 

In deriving these equations, the boundary condition was 
used; it states that 1(0) = /c(0) as no scattering has 
been incurred yet at the interface. According to energy 
conservation, the incoherent intensity Ijj (diffusive) is 
thus 

I D {x) = I(x) - I c {x). (13) 

When there is no absorption, the decay constant a 
is expected to vanish and the total intensity will then be 
constant along the propagation path. Then, the coherent 
energy gradually decreases due to random scattering and 
transforms to the diffusive energy, while the sum of the 
two forms of energy remains a constant. This scenario, 
however, changes when localization occurs. Even without 
absorption, the total intensity can be localized near the 
interface due to multiple scattering. When this happens, 
a does not vanish. The transport of the total intensity 
may be still described by Eq. I|10|l . and the inverse of a 
would then refer to the localization length. 

2. Diffusion theory 



side of a sample and measuring the transmission across 
the sample on the other side, may be influenced rather 
significantly by other non-localization effects such as re- 
flection, in the present paper we use an alternative setup, 
in line with the previous studies |32j, |4fJ : Put the source 
in the middle of a scattering sample which takes a circular 
shape. In this case, a solution to the diffusive intensity 
is given as [32| 

Id = -\\n{r/R) + I b , (17) 

■Kit 

where R is the samole radius, i. e. the size of the cluster 
of the scatterers, and lb is the intensity at the boundary. 
In practice, l t can be regarded as a free parameter to 
fit the numerical simulation. When the diffusive energy 
is dominant, the total energy would be expected to also 
follow the behavior characterized by Eq. (|17fl . 

III. THE RESULTS AND DISCUSSION 

A conceptual set up of the system for discussion is pic- 
tured by Fig. Here it is shown that a source is located 
at the center of a sample circled by the broken line. The 
small circles inside refer to the dipoles as the scatterers. 
The receiver is placed at various spatial positions along 
the radial direction to record the transmitted waves. 



As discussed above, as waves propagate, the coherent 
energy gradually decreases and the incoherent (the so 
called diffusive waves) progressively grows. In the ab- 
sence of localization, the diffusive intensity (or energy) 
will eventually surpasses the coherent energy at a cer- 
tain point, then becomes the dominating component in 
the wave transmission. This intuitive physical picture 
has been nicely phrased in the very excellent textbook 
by Ishimaru |39| . 

The diffusive waves are expected to follow the diffusion 
equation which can be described by 



(dv 2 ) i u = -mr- F, 



(14) 



where D is the diffusive constant and S denotes a 
strength factor which is related to the total emitted en- 
ergy (Jo) and transport mean free path (It}- In the steady 
situation, this equation is simplified as [39| 

V 2 I D = -^6(r-r s ). (15) 

The boundary condition for the diffusive waves is |39| 



, , 2l t dl 
3 on 



= 0, 



(16) 



with n being the outward normal at the boundary. 

As discussed in [4(| , the traditional wave of measuring 
for localization effects, that is putting the source at one 
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FIG. 1: Conceptual layout of the system and the simulation. 



A. Parameters in the simulation 

Unless otherwise noted, the following parameters are 
used in the numerical simulation: the non-dimensional 
damping rate, bo/ujo = 0.001 and the interaction cou- 
pling, b = 0.001. Without notification, the filling factor 
(j3) is taken as 6.25; the filling factor is defined as the 
number of dipoles per unit area. The number of ran- 
dom configurations for averaging is taken in such a way 
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that the convergency is assured. In the calculation, we 
scale all lengths by a length I such that k l = 1, and 
frequency by u)q. In this way, the frequency always en- 
ters as k/ko. We find that all the results shown below 
are only dependent on parameters b, bo/coo, and the ra- 
tio uj/luq or equivalently k/ko. Such a simple scaling 
property may facilitate designing experiments. In the 
numerical computation, we take c = 1 for convenience. 
The total normalized wave at a spatial point is scaled as 
T(r) = E(f) / E (r), with E being the direct wave from 
the source, so that the trivial geometric spreading effect 
is naturally removed. When compared with diffusion the- 
ory, the diffusive energy would have to be multiplied by 
the spreading factor, and then is Idt. 



B. Numerical results and discussion 
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FIG. 2: The transmission versus the frequency. 

First we plot the wave transmission versus frequency 
to locate the strong localization region. The result is 
shown by Fig-EI Here totally 1964 dipoles are considered, 
amounting to a sample radius of about 10. The receiver 
is located outside at about 1.2 radius of the sample. We 
see that there is a very strong inhibition region ranging 
roughly from 1.001 up to 1.012. Outside this regime, 
the inhibition is not obvious. We have tried our current 
computing facility to its extremes with the largest possi- 
ble number of dipoles, but we still failed to see obvious 
inhibition outside the regime. The width of the strong 
inhibition window mainly depends on the the filling fac- 
tor. We choose five frequencies, indicated on Fig. |3 to 
study the localization or diffusion behaviors of the wave 
intensity. 

In Fig. |21 we plot the transmission behavior of vari- 
ous intensities as a function of wave travelling distance. 
Enough average over the random configurations of the 
scatterers is taken to ensure the stability of the results. 
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FIG. 3: Transport behavior: (a) the total, coherent and dif- 
fusive intensity as a function of travelling distance; (b) the 
comparison of the numerical results of the total energy with 
that from the diffusion theory; (c) the comparison of the nu- 
merical results of the diffusive energy with that from the dif- 
fusion theory. We considered 1964 dipoles to form a sample 
whose radius is roughly 10. 



The maximum number of the average is 1200. The obser- 
vation can be summarized as follows. (1) It is clear from 
the left panel that within the strong inhibition window, 
the transmission follows perfectly the expected localiza- 
tion behavior; the total and coherent intensity decays ex- 
ponentially with slight different slopes respectively, while 
the diffusive intensity increases initially, then saturates, 
followed by a decay afterwards. These are shown by (al). 
(2) As the frequency moves towards the edge of the strong 
inhibition window, the behavior starts to change more 
and more severely. As the frequency increases gradually, 
the total energy starts to behave more like the diffusive 
waves, although the diffusive waves are not yet dominant 
(such as in (a2)). But overall speaking, the transport be- 
havior of the total energy behaves more like the diffusive 
waves when the latter are prominent, as expected and 
shown by (a3) and (a4) . (4) At frequencies for which the 
diffusive intensity is dominant, the behavior of the diffu- 
sive energy is qualitatively in agreement with the predic- 
tion of the diffusion theory, such as in the cases described 
by (c3), but in certain cases the agreement is good both 
qualitatively and quantitatively as seen in (c4). Over- 
all speaking, however, the agreement is less encouraging 
compared with the simulation for another 2D system il- 
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lustrated in [S^]. And the rapid fluctuation of the diffu- 
sive energy shown in [32l is absent here. Another note is 
worth making here. In the authors also compared 
the results within a strong localization region. The expo- 
nential decay revealed in their figures, for example Fig. 9, 
is only possible when the cylindrical geometrical spread- 
ing factor 1/r is taken out. It seems to us that this factor 
has not taken into account in their discussion. (5) Just 
by looking, we are tempted to conclude that there is in- 
deed a fundamental difference in the transport behaviors 
between within and outside the strong inhibition region. 
This conclusion accords qualitatively with j^]. (6) We 
also observe that outside the strong inhibition regime, 
the coherent intensity dose not necessarily behave as an 
exponential decay, exemplified by (a3) and (a4). This is 
in contrast to usual expectation. A plausible cause is the 
finite size effects. In this situation, the mean free path is 
hard to estimate, unlike in other systems p9ll40|. 
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FIG. 4: Dependence on the sample size of the transport be- 
havior. The following parameters are taken. (1). k/kO = 
1.013; (2). Filling factor = 6.25; (3). The number of average: 
400; (4). Fig(al to cl) : dipole number: 491, Sample ra- 
dius: 5; Fig(a2 to c2): dipole number: 982, Sample radius: 7; 
Fig(a3 to c3) : dipole number: 1473, Sample radius: 9; Fig(a4 
to c4) : dipole number: 1964, Sample radius: 10; Fig(a5 to 
c5) : dipole number: 2455, Sample radius: 11. 

From Fig. |31 we see that at around k/k = 1.013, just 
outside the upper edge of the strong localization bound- 
ary, the diffusive waves are dominant and the transport 
behavior tends to follow the prediction from the diffusion 
theory. We would like now to explore the robustness of 



such an agreement. The results are shown in Fig.^J Here 
it is shown that as the size of the sample is increased, the 
diffusive intensity become more and more prominent and 
meanwhile the agreement between the diffusion and exact 
numerical results is improved, as shown by the diagrams 
from the top downwards. When we consider another fre- 
quency fc/fco — 1-014, the agreement seems not improving 
as the sample size is enlarged, although the two frequen- 
cies are close. The comparison is shown in Fig. [5J The 
best agreement is when the sample size is 7, referring to 
(b2) and (c2). From these results, one may say that the 
agreement between the diffusion theory and the exact nu- 
merical results may not be robust, at least in the present 
system. We have to note that the results from the dif- 
fusion theory are approximate. The deviation from the 
diffusion theory is therefore not out of expectation. 
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FIG. 5: Dependence on the sample size of the transport be- 
havior. The following parameters are taken. (1). k/kO = 
1.014; (2). Filling factor = 6.25; (3). The number of average: 
400; (4). Fig(al to cl) : dipole number: 491, Sample ra- 
dius: 5; Fig(a2 to c2): dipole number: 982, Sample radius: 7; 
Fig(a3 to c3) : dipole number: 1473, Sample radius: 9; Fig(a4 
to c4) : dipole number: 1964, Sample radius: 10; Fig(a5 to 
c5) : dipole number: 2455, Sample radius: 11. 

We also considered the case for another filling factor at 
3.00. The general features are agreeable with the above 
case with 6.26. For brevity, we just show partial results. 
Fig. [S] shows the transmission versus the frequency. The 
parameters are shown in the figure caption. Compared to 
the case with filling factor of 6.25, we see that decreasing 
filling factor only tends to reduce the strong inhibition 
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regime while all other features remain unchanged. 



Filling Factor = 3.00 
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k'ko 

FIG. 6: The transmission versus the frequency. The param- 
eters are as follows. (1) 2150 dipoles in a circular sample, 
amounting to a sample radius of 15; (2) the filling factor = 
3.00; (3) the receiver is positioned at 1.2 radius of the sample. 

Again we choose five frequencies to inspect the trans- 
port behavior, and show the results in Fig. [7| The gen- 
eral conclusions drawn from Fig. [3] remain. However, 
there are differences. For the frequency 1.004 which is 
within the strong localization regime, the transport be- 
havior fully complies with the expectation. The diffusive 
intensity increases, saturates, and then decreases to take 
over as the dominating part in the total transmission. We 
re-plot fig. [3 (al) in Fig. 8 to show this behavior. The 
cross over point when the diffusive part becomes domi- 
nating is represented by the vertical line in the figure. 

Lastly, we have computed the fluctuation of the nor- 
malized transmission intensity versus frequency. The re- 
sults are shown in Fig.|5J The parameters are as follows: 
filling factor (Number of dipole/area) = 6.25; two sam- 
ple sizes are 5 and 10 respectively; the receiver is at 1.2 
times the radius; the number of average is 100. The re- 
sults indicate that the fluctuation is nearly zero with the 
localization regions. Outside these regimes, the fluctua- 
tion is strong. The results also show that the localization 
regime is the same for different sample sizes. These fea- 
tures hint that there may have a localization transition 
in 2D random media, in accordance with the previous 
results [Hill- 

IV. SUMMARY 

In this article, the transport properties of electromag- 
netic waves in a two-dimensional random system. Some 
general properties of the transport phenomena are elab- 
orated. For certain ranges of frequencies, strongly local- 
ized electromagnetic waves have been observed in such a 
simple but realistic disordered system. The spatial be- 
havior of the total, coherent and diffusive waves is ex- 
plored within, near the edge of and outside the localiza- 
tion regime, and are investigated for different parameters 
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FIG. 7: Transport behavior: (a) the total, coherent and dif- 
fusive intensity as a function of travelling distance; (b) the 
comparison of the numerical results of the total energy with 
that from the diffusion theory; (c) the comparison of the nu- 
merical results of the diffusive energy with that from the dif- 
fusion theory. We considered 2150 dipoles to form a sample 
whose radius is roughly 15. 
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FIG. 8: Re-plotted results in Fig.^al). 



such as filling factor and system size. The results show 
that within the localization regime, waves are trapped 
near the transmitting source. Meanwhile, the diffusive 
waves follow an intuitive but expected picture. For the 
cases that the frequencies are near the boundary of or 
outside the localization regime, the results of diffusive 
waves are compared with the diffusion approximation, 
showing less encouraging agreement as in other systems 
[32j |. Furthermore, the results tend to suggest that the 
study of the fluctuation behavior in the transmission may 
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also help identify the localization regime. 
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